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ing finite energy per unit lengtli is constructed numerically in an extended Abelian 



, A new class of twisted, current carrying, stationary, straight string solutions hav- 

(N 

' Higgs model with global SU(2) symmetry. The new solutions correspond to defor 

Mh ; 

, mations of the embedded Abrikosov-Nielsen-Olesen (ANO) vortices by a twist - a 
relative coordinate dependent phase between the two Higgs fields. The twist induces 

' a global current flowing through the string, and the deformed solutions bifurcate 

■ with the ANO vortices in the limit of vanishing current. For each value of the wind- 



ing number n = 1,2... (determining the magnetic flux through the plane orthogonal 



(N 
O 

f-*) I to the string) there are n distinct, two-parametric families of solutions. One of the 



continuously varying parameters is the twist, or the corresponding current, the other 

Oh. one can be chosen to be the momentum of the string. For fixed values of the mo- 

CD ^ 



mentum and twist, the n distinct solutions have different energies and can be viewed 
as a lowest energy "fundamental" string and its n — 1 "excitations" characterized 



^ . by different values of their "polarization". The latter is defined as the ratio of the 

angular momentum of the vortex and its momentum. In their rest frame the twisted 
vortices have lower energy than the embedded ANO vortices and could be of con- 
siderable importance in various physical systems (from condensed matter to cosmic 
strings) . 
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I. INTRODUCTION 

It was almost 50 years ago that Abrikosov predicted the existence of vortex hnes Q| in 
type II superconductors in the framework of the Ginzburg-Landau (GL) theory, which has 
then been confirmed many times by experiments in various superconducting materials. The 
importance of such vortex lines or strings in high energy physics has been recognized by 
Nielsen and Olesen, who have shown the existence of classical solutions describing straight, 
infinitely long strings with finite energy per unit length in a relativistic field theory - in 
the Abelian Higgs model j^. Ever since their discovery, vortices have been the subject of 
intense studies, since they find numerous applications in many domains of physics including 
various branches of condensed matter theory (such as superconductivity, superfluid ^He sl), 
high energy physics and cosmology. In a more general context vortex lines belong to the 
class of objects known as "topological defects". Topological defects (monopoles, vortices, 
domain walls, etc.) result from the spontaneous breaking of global or local symmetries of the 
underlying field theory; we refer to some of the excellent reviews |4|. There is an interesting 
class of models in which both global and local symmetries are simultaneously broken, in 
the most economical way, with the minimal number of scalar fields. These models, dubbed 
semilocal models, have some remarkable features; for a detailed review see Ref. ^. Perhaps 
the most interesting point is that semilocal models exhibit stable vortices, despite the first 
homotopy group of the vacuum manifold being trivial [6,Q]. 

A particularly interesting and simple class of such semilocal models is the extended 
Abelian Higgs model (EAH), which contains electromagnetism coupled to N charged com- 
plex scalars transforming as a vector in the fundamental representation of the global SU(N) 
symmetry group. The simplest nontrivial EAH model is the SU(2) symmetric (N=2) one, 
which case is made quite interesting by noting that it represents the bosonic sector of the 
Electroweak theory in the limit where the Weinberg angle, = 7r/2 A non- relativistic 
version of this model is also used in condensed matter theory, in the description of un- 
conventional superconductors for example P|. The Abrikosov- Nielsen-Olesen (AND) vortex 
solutions P], characterized by an integer winding number, n, determining their magnetic 
flux, and by the mass ratio, = ms/m^ (mg resp. denoting the mass of the scalar 
resp. vector fields), can be embedded into this theory. The parameter (3 also distinguishes 
superconductors of type I (/5 < 1) from type II (/3 > 1). The case /5 = 1 is quite special, for 
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here, instead of there being a unique vortex solution satisfying the first order Bogomol'nyi 
equations , there is a continuous family of them having the same energy [3| . The stability 
of the ANO vortices embedded into the EAH models was examined in Ref. 0, and it was 
shown there that they are stable only ii (3 < 1. For (3 > 1 the embedded ANO vortices 



have been found to ex 
spreading" instability 



libit a single unstable mode that can be interpreted as "magnetic 
3, llll • The Standard Model counterparts of the embedded ANO 
solutions are the Z-strings, which have been shown to be stable for sin^6'-w ^ 0.9 jl^ . 

In this paper we give a detailed description of a new class of twisted, straight string 
solutions in the simplest nontrivial SU(2) symmetric semilocal EAH model, some of our main 
results having been announced in Ref. 1^. The new solutions are translationally symmetric 
in the z = direction. The main feature distinguishing them from the previously known 
string or equivalently vortex (when the string is viewed in the plane) solutions is 

due to the presence of a relative phase (twist) between the two scalars, depending linearly 
on z. The twist induces a persistent global current flowing along the z-direction. It turns 
out that there is a maximal value of the twist, depending on j3. For this maximal twist the 
current vanishes, and the twisted solutions bifurcate with the embedded ANO vortex. In 
fact the above mentioned "magnetic spreading" instability of the embedded ANO vortex 
in the EAH model can be viewed as a signal of bifurcation with a new branch of current 
carrying solutions. 

For any fixed value of (3 > 1, and for a given winding number, n, the static twisted vortices 
form n distinct one parameter families, labeled by the twist (or equivalently by the current) 
and by an integer u = 1, . . .n. For a fixed value of the twist, the n different solutions have 
different energies and can be viewed as a lowest energy "fundamental" string and its n — 1 
"excitations". It is quite interesting that the energy of the fundamental twisted solutions is 
lower than that of the corresponding ANO ones. They could therefore have better stability 
properties, and in fact we expect the fundamental n = 1 twisted string to be stable in the 
N= 2 EAH model for any value of (3 > 1. 

In the simplest case the twisted vortices are static. But the relative phase of the two scalar 
fields may also have a linear dependence on the time coordinate, t = x^, in which case the 
vortices become stationary. These solutions can be obtained from the static ones by Lorentz 
boosts in the t, z plane. The stationary strings are surrounded by an exponentially screened 
radial electric field, and they also possess a nonzero momentum, as well as a nonzero angular 
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momentum oriented in the 2;-direction. For any fixed value of /3 > 1, and for a given winding 
number, n, the twisted, stationary vortices form n distinct two-parameter famihes which 
can be labeled by their momentum and current. The discrete parameter, u = 1, . . . n, which 
distinguishes between the n different families can be interpreted as a kind of "polarization" , 
since it determines the ratio of the momentum and the angular momentum of the string. 
Due to the Lorentz contraction, the energy of the stationary strings is higher than that of 
the static ones. 

Our solutions can also be viewed as superconducting cosmic strings in a "minimal" model, 
with a single U(l) gauge field, unlike in the prototype U(l)xU(l) model of superconducting 
strings introduced by Witten Q]. In contradistinction to the superconducting strings in 
Witten's model, the twisted semilocal strings are strongly localized in the transverse di- 
rection, which guarantees the finiteness of their energy. In addition, their superconducting 
current is global and not local, and it can be arbitrarily large, unlike the situation in Witten's 
model [15^ . 

Most previously known semilocal string solutions are static and untwisted, with the no- 
table exception of the solutions of Abraham [16], which have however infinite energy per 
unit length (for n = 1). We should also like to mention that in Ref. in a theory similar 
to the EAH model, straight, twisted string solutions have been obtained in the description 
of neutral two-component plasmas. Finally it should be pointed out that the non-abelian 
counterparts of our twisted strings should exist in the Electroweak model in the form of 
"twisted Z-strings" in analogy to the ordinary Z-strings. 

The plan of the paper is the following: In Section II we introduce the EAH models. In 
Section III the stationary, z-translationally symmetric Ansatz is given, and the EAH theory 
is dimensionally reduced to two dimensions corresponding to the plane orthogonal to the 
symmetry axis. Section IV is devoted to the classification of the stationary and twisted string 
configurations according to their properties under Lorentz transformations. In Section V 
the theory is further simplified by considering rotationally symmetric configurations in the 
plane. We discuss the field equations, describe the local behaviour of their solutions in the 
neighborhood of the singular points. Section VI contains a short review of some relevant 
known vortex solutions: the classical ANO vortices, the "skyrmion" solutions satisfying the 
Bogomol'nyi equations, and the chiral, charged vortices of Ref. Section VII is the main 
part of the paper, it contains a detailed description of the new solutions. Some concluding 
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remarks are given in Section VIII. 

II. EXTENDED ABELIAN HIGGS MODEL 

The extended Abelian Higgs (EAH) model considered below is a theory of N complex 
scalars, = (0i,...0n), with an SU(N) global symmetry and with their overall phase 
gauged. Such a theory possesses an SU(N)giobaixU(l)iocai symmetry algebra by construction. 
In this paper we shall concentrate on the simplest case corresponding to N=2, mostly because 
of its relationship to the electroweak model. The 4 dimensional (4D) action of the theory is 
given by 

S = F^^F'^'' + {D^^yD^'<^ - ^ ($t$ - r]^^ (fx . (2.1) 

Here F^^, = d^A^ — dyAi^, is the gauge field tensor, D^(^ = dfj,^ — igA^^ is the gauge covariant 
derivative, g denotes the coupling (charge), and the signature of the flat Minkowskian metric 
used here is (+, — , — , — ). Due to the form of the potential the scalar field has a vacuum 
expectation value, t], and the U(l) gauge symmetry is spontaneously broken. The physical 
spectrum consists of a vector particle with mass rriy = V2gr], of 2(N - 1) Nambu-Goldstone 
bosons, and of a Higgs scalar with mass rris = \/2Xrj. 

It is convenient to scale the fields and coordinates as: $ 1]^, — > rjA^, —>■ / {grj), 
so the rescaled fields and coordinates are dimensionless. In terms of the rescaled quantities 
the action ()2.1|) becomes: 

5=^1 Cd'x=^^ jd'x |-iF^,F'^'^ + (D^$)tD^$-^(<|.t$-l)2| , (2.2) 

where now D^^ = 9^$ - 2^^$, and /3 = X/g"^. 

The fields transform under the U(l) gauge symmetry as 

A^^ A, + d^K{x), $^e^A(^)<l>, (2.3) 

while the global SU(N) symmetry acts on the fields in the following way: 

A,^A^, $ ^ e^C'"^")*, (2.4) 

where (a = 1, . . . , N) denote the generators of the SU(N) Lie algebra in the fundamental 
representation (for N=2 we choose to be the Pauli matrices), and the da are constant 
parameters. The corresponding Noether currents are 

jI = -^{(D^$)tT"$ - $tT"D^$}, (2.5) 
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with a = (0, a) and = {1,t°') being generators of SU(N)xU(l). The conservation of the 
currents ()2.5p follows from the Euler-Lagrange equations for the action ()2.H) . which can be 
written as: 

d'Fp, = -jI = - , (2.6a) 

DpDP^ = -1)^ , (2.6b) 

where |$p = 

We shall also need the energy momentum tensor induced by varying the metric tensor: 

dC 

T'. = 29'''g^-5'.^- (2.7) 
With C given by ()2.2|) this evaluates to 

T'; = + (D^<l>)tD^$ + {D^<t>yD^<^ - d^'^C . (2.8) 

The energy density reads 

n = \Do^' + 1 A$r +^{m'-^Y + \E'+ , (2.9) 

where the Latin indices take the values = 1,2,3, |-Di$P = {Di^)'^ Di^^ (summation 

over the index i is implied), the electric and magnetic fields, E and 5, are defined as Ei = Foi 
and Bk = eijkFij/2, respectively. 

III. DIMENSIONAL REDUCTION TO 2 DIMENSIONS 

Modulo gauge and global symmetry transformations the ground state of the theory (j2.2p 
(i.e. the solution of zero energy) is given by = 0, $ = $0) where $o is a constant 
complex N-vector, (|$oP = !)• Thus for the case considered in this paper, N = 2, the 
vacuum manifold is a 3- sphere: 

V = {$ G I $t$ = 1} = . (3.1) 

S symmetric, finite energy solutions of the EAH 
. These strings are infinitely long and their total 
energy is infinite, nevertheless they can be of physical interest if their energy per unit length 
is finite. In what follows we shall be considering such strings in the simplest case when 



As is well known, there are no spherical 
model, but it admits stringlike solutions 
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they are stationary and translationally invariant along the z = axis. In this case the 
most general Ansatz for the vector field and for the two (upper and lower) complex 
components 0i and 02 of the Higgs field $ can be written as 

= {Ao{xi,X2) ,A3{xi,X2) ,Ai{xi,X2)) = {Ao,,Ai), a = 0,3, i = l,2, 
01 = X2), 02 = f2{xi, x2)e*(-»*+-^^) , (3.2) 



Ua, are real parameters. This Ansatz 
applied to the case of two commuting 



where /i,/2 are complex functions and (cjq, i-t^a)^ 
easily follows from the general results of Ref. 
symmetry transformations. The most important feature of the Ansatz ()3.2j) for the present 
paper is that there is a nontrivial, coordinate dependent relative phase between the two 
components of the SU(2) doublet 0i and 02- The change in the phase of the scalar fields 
induced by infinitesimal space-time translations in the t, z directions is compensated by 
gauge transformations (followed by global SU(2) rotations). In other words a space-time 
translation moves the fields along gauge orbits. The relative phases, Ua, can be interpreted 
as a relative rotation, uq, resp. twist along the z-axis, u^, between the two components of 
the scalar field (0i,02). One should also note that the Ansatz ()3.2|1 breaks the originally 
present global SU(2) symmetry to U(l). 

The vacuum manifold ()3.1|) is a three sphere, which is simply connected (7ri(S'^) = 1), 



3, 



because of the 



and so there are no topological string solutions. However, as explained in 
spontaneous breaking of the U(l) gauge symmetry the space of finite energy configurations 
is not labeled by the vacuum manifold V, but rather by the gauge orbits, which space is 
not simply connected: vri(U(l)iocai/l) = ^- This implies that configurations with different 
U(l)iocai winding numbers are separated by infinitely high energy barriers. 

For field configurations given by Ansatz ()3.2|) the action ()2.2|) dimensionally reduces from 
4D to 2D, 5 = / Sdx^dx^ where 

S = V' J rf'x{A„A°|0j2+^-(^^-2A,)|02|2-|A0ar + ^5.Aaa,A°-i4-| (I0J2 _ 1)2} , 

(3.3) 

where a = 1,2; the index a = (0, 3) is raised by the residual Minkowskian metric, i.e. 
A°'Aa = Aq — A"^, while the 2D metric of the reduced theory in the {x^,x^) plane is Euclidean; 
summation over the indices a,i,j = 1,2 is understood. 
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The equations of motion of the reduced theory ()3.3|) in the (xi, X2) plane take the form: 
AA^ = 2AM^ - ^^M' , (3.4a) 



diFij = t{<PaD,(j)^ - <p,Dj<Pa) , (3.4b) 

AA0a = - l)0a " ^a^>a " C^"(c^a " 2A„)025a • (3.4c) 

For configurations satisfying the above Gauss law- type equations ()3.4b ). the energy per unit 
length of the string is given by £^ = rfE, where the " dimensionless" energy per unit length, 
used throughout this paper can be expressed as: 

E=j d'xT\ = j d^x[uj^{uj^ - + + \D4a? + f (|0aP - l)'} , (3.5) 

where B stands for the z-component of the magnetic field (orthogonal to the (xi, X2) plane); 
B = B3 = F12 = eijdiAj, and summation over the repeated a = 0, 3-indices (not raised with 
the residual Minkowskian metric) is understood. 

The Noether current associated to the residual U(l) symmetry of the ansatz ()3.2|) actually 
coincides with the 3-rd isotopic component, j^, of the global SU(2) current in Eq. ()2.5p . 
Introducing ^"''^J^ = — i(0aD^0^ — ^^i^^^a), where now no summation is implied over a = 1, 2, 
one has = ^^^Jfi — On the other hand, the local U(l) current in Eq. (j2.5p can be 

written as j° = ^^^J^ + '-^-'J^. Since, by virtue of the field equations ()2.6|) . for fields decaying 
sufficiently fast at infinity the integral of j° over the two-plane vanishes, it follows that the 
integral of the current can be expressed as 

j Sx]l = 2 j d^x = "2 / ^^^-^^ •= 
where the "string worldsheet current", X^, is 

Ia= d'^X Aa(j)l(f)l = / d'^x{uJa - Aa)(j)24>2, (3.7) 



whose components are the conserved Noether charge per unit length of the string, Xq, and 
the total current fiowing along the string, X3. 

Due to the translation symmetry of the Ansatz ()3.2j) there is a conserved momentum, P, 
which is expressed (per unit length) as 

P = I d^xT^^ = 2U0I3, (3.8) 
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and for configurations which possess an additional rotational symmetry in the plane (which 
is the case for all solutions presented in this paper) one also has a conserved angular mo- 
mentum, J (per unit length), given by: 

J = j (i^xT^. (3.9) 

Note that just as for the energy per unit length, £ = rfE, to obtain the dimensionful 
momentum and angular momentum one should multiply P and J hy rf . 



A. Energy bounds and Bogomol'nyi equations 

Let us recall briefly the method, established by Bogomol'nyi [10], to find a lower bound 
for the energy functional We first introduce the total magnetic flux : 

j (fxB = j Ajdx^ = 27m , (3.10) 

with n being an integer winding number. The flux ()3.10|) is quantized and topologically 
conserved, due to the non-trivial topology of the gauge orbits. Thus configurations with 
different values of n are separated by infinite energy barriers even though the potential 
energy cost can be arbitrarily small, since the Higgs field is topologically trivial 6]. In 
what follows, we shall refer to n as the winding number or number of magnetic flux quanta. 
Without any loss of generality we can assume that n > 0. Next, using the identity 

^|A$ + = + e,,.9,($tD,$) - 5|<l>|2 (3.11) 

and the expression dSHI) for Ja, the energy functional ()3.5|1 becomes 

E = 27m+\u;M+j (fx\^l\D,<i>+ie,,Dj^^ + l{B+m^-lf + l{p-l){\<^\^~lf^ . (3.12) 

It can now be seen that for the special value (3 = 1, the energy is minimized to 2?™ + \uJaTa\ 
when the following first order equations, the Bogomol'nyi equations, are satisfied: 

Di^ + ieijDj^ = 0, 
B + |<l>|2 -1 = 0. (3.13) 

If in a fixed sector of the phase space labeled by the winding number n, the value of \u!aXa\ 
can be maintained fixed when varying the action ()3.3|) . then solutions of the Bogomol'nyi 
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equations ()3.13|) will solve the variational Euler-Lagrange equations ()3.4|) . In the static and 
untwisted case, when Ua = 0, this is automatic, but for Ua ^ this is not true in general, 
except for solutions with UaUj"' = discussed in the next Section. 

The known static and untwisted {ua = 0) semilocal solutions are either the embedded 
ANO vortices, with Ai = A(ANO)i and $ = (f)(ANO)^o, where $o is a constant SU(2) doublet 
of unit norm, $o$o = 1, or the so-called "skyrmion" 0, 19 1 solutions. The latter fulfill the 
Bogomol'nyi equations ()3.13p . and for any winding number, n, there exists a 4?T,-parameter 
family of such solutions describing multi-vortices. 2n of these parameters can be identified 
with the positions of the vortices (n points in the plane), and every vortex located at a point 
p is labeled by a complex parameter, Wp, whose modulus determines the width of the flux 
tube at point p, and whose argument corresponds to the orientation of the vortex in the 
"internal" space. 



IV. LORENTZ CLASSES 

We notice at this point that our coordinate system is not yet completely fixed, since one 
can perform boosts along the 2;-axis: 

= (t, z, xi, X2) x'^ = {i, z, xi, X2), (4.1) 

where 

t = tcosh(7) — 2sinh(7), 2; = icosh(7) — t sinh(7), (4-2) 
with 7 being the parameter of the boost. The fields then transform as 

A^{x)dx'' = i^(x)rfx^, $(x) = 4>(x) , (4.3) 

and the form of our Ansatz ()3.2j) will be preserved upon replacing t, z hj t,z and {ujq, U3) 
by 

ujq = LUo cosh(7) — U3 sinh(7) , uj^ = 003 cosh(7) — uq sinh(7) . (4.4) 

We thus see that Lorentz boosts in the (t, z) plane preserve the structure of the ansatz ()3.2p . 
and that the Ua transform as components of a Lorentz (co)vector. This means that all pairs 
of values of cUq belonging to the same orbit of the Lorentz group (see Fig. ^ are equivalent. 
As a result, instead of two independent parameters {ujo, uj^) we have in fact only one, which 
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can be chosen to be the Lorentz-invariant combination 



= —UJ°'Ua 



^3 - ^0 • 



(4.5) 



Therefore the space of solutions decomposes into three classes that can be labeled by the 



possible Lorentz types of the length of u"^, as proposed by Carter 

= null or chiral case 
< time-like or electric case 
> space-like or magnetic case 



(4.6) 
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FIG. 1: Phase space of the stationary and translationally symmetric strings. Solutions correspond- 
ing to a Lorentz orbit with constant u'^ = w'^ — loq are related to each other by the boost symmetry. 
The new superconducting strings correspond to the u)^ > regions. 

These three classes also correspond to the Lorentz types of the string worldsheet current 
la in (!3.7j) . We now discuss some relevant properties of Eqs. (j3.4p corresponding to the 
three Lorentz classes. 

1. The chiral class, UaUj"" = 0. All known semilocal solutions obtained up to now belong 
to this class. From Eqs. ()3.4b ) it follows that A(^^) := u'^Aa satisfies 

AA(,) - 2A(^)0„0, = 0. (4.7) 
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Now the standard positivity argument for Eq. ()4.7p . 

J rf\A(^)AA(^) = - J d^xid^A^^)^ = 2 J rf2a;(A(^))Va0a, (4-8) 

implies that the only regular and bounded solution of Eq. ()4.7|1 in the 2-plane is the 
trivial one, = 0, i.e. ujqAq = u^A^. Therefore the two Gauss-law type equations 
()H.4k ) reduce to a single one. This also implies that the worldsheet current is isotropic, 
Xq,X" = 0. Furthermore, since A^^A^ = 0, Eqs. ()3.4b .c) reduce to the equations of 
motion of the untwisted and static semilocal vortices, and they decouple from the 
remaining Gauss law. Therefore, to any static, untwisted vortex solution {Ai,(f)a) of 
Eqs. flH.4b .c) there corresponds a twisted, stationary chiral vortex obtained by solving 
in addition the Gauss-law equation ()3.4k ). It is important to note that the solutions 
Aa of the Gauss-law equation ()H.4b .) do not always fall off sufficiently fast at infinity. 
As a result, not all of the stationary chiral solutions will have finite energy. 

As already mentioned, the known static and untwisted semilocal solutions of Eqs. 
()H.4b .c) are either the embedded ANO vortices, or the "skyrmions". Bounded solutions 
of the related Gauss-law equation ()3.4b ) exist only for the "skyrmions", but their 
energy for winding number n = 1 is nevertheless divergent due to the slow falloff of 
the fields A^ (lack of localization). This is an alternative way to rederive the main 
result of Ref. where stationary and charged vortices satisfying the Bogomol'nyi 
equations were found, to which we shall refer to as "spinning skyrmions" . 

2. The space-like (or magnetic) class, UaUj" < 0. In fact all new solutions to be presented 
in this paper fall into this class. Now contracting Eqs. ()3.4b ) with Ua where u'^Ua = 
one finds that A(^) := u°'Aa = uj^Aq — uqA^ satisfies 

AA(^)-2A(^)0j, = O. (4.9) 

The same positivity argument as above implies again 74(^) = 0, i.e. the electric po- 
tential, Aq, is completely determined by the component of the magnetic one, A^. 
Exploiting the boost-symmetry of the Ansatz (j3.2|) . one can always achieve in this 
case by a choice of the reference frame that ujq = and Aq = 0. We note that the 
worldsheet current is spacehke, X^^^" < 0. 

3. The timelike (or electric) class, 00^00°" > 0. In this case there is such a reference frame 
where cus = and ^3 = 0. The worldsheet current is timelike, Xq,X° > 0. It is unlikely 
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that finite energy solutions exist in this class, since from Eqs. ()3.4b ) we obtain the 
asjTiiptotic (p^ = xl + xl oo) equation A/2 = —uJaUj"^f2 (assuming that fi — > 1, 
~^ 0), which implies that the function I/2I tends too slowly to zero to make the 
energy ()3.5|) convergent. 

V. THE CYLINDRICALLY SYMMETRIC EQUATIONS 

A. Field equations 

In this section we shall examine solutions with cylindrical symmetry, for which we may 
choose a gauge such that we are lead to the following Ansatz in polar coordinates: 

Aq = ujottoip) ,Ap = 0,A^ = na{p) , A3 = u;3«3(p) , 

01 = /i(p)e'"^ 02 = /2(p)e*'"'^e*('^»*+-^^) , (5.1) 

where m = 0, . . . n — 1. Furthermore, as shown in the previous section, the electric potential 
is either given by Aq = A3 in the chiral case {\ujq\ = Ic^sl), or by Aq = uqA^/us in the 
magnetic case, i.e. in both cases one can take ao(p) = ct3{p). The energy functional assumes 
the form: 

/2 . n 



E =2n r pdpT\ = 2n H pdp\Un'^ + H'a'A +n'^-^-^f^ 
Jo Jo I ^ ^ P P 

+ ^"'p^' ii + f[' + r.' + + m - a,f] + § (1 - i/p)^} , 



(5.2) 



where = u^ + ojI, |/p = /i +/|. A Bogomol'nyi-type rearrangement of the energy (|5.2p 
yields the following: 

00 

E = 27rn + \uj^\Q + 7r(/3 - I) j pdp (l - \f\^Y+ 

(5.3) 





00 







where the quantity Q, 



00 00 

2 _ o„ / „J„ „ /2 



Q = 2ttJ pdp{l- as)f^ = 2nJ pdp a^ff , (5.4) 
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determines the vortex worldsheet current 



la = U)aQ, (5.5) 

and in terms of which the momentum ()3.8|) and the angular momentum can be ex- 
pressed as 

P = 2U0U3Q , 

J = —2ujqvQ , where v := n — m = 1, . . . ,n . (5-6) 

The parameter Q is Lorentz- invariant, since it enters the invariant relation ()5.5p between 
Lorentz-vectors. On the other hand, the quantities Xo,P, J are not Lorentz- invariant. In 
particular, in the spacelike case, if uj°'uja < 0, there exists a rest frame where ujq = 0, 
therefore in this frame Xq = P = J = 0. Using Eq. ()5.6|) . we can define in any nonstatic 
frame a "helicity" or "polarization" given by the ratio 

^ = --. (5.7) 

P UJ3 

which can be extended as a limit to the static reference frame as well. From now on we shall 
refer to u as helicity number. In the rest frame, where uj^ = uj, we define the "rest frame 
current" as 

X:=ujQ. (5.8) 

The value of the energy in the rest frame is given by Eq. ()5.3|) with |u;p = a;^. 
Eq. ()5.3|) implies the following inequality for (3 > 1: 

E>27m+ {u^ + ul)Q . (5.9) 

It is worth noting that the existence of this bound does not automatically guarantee the 
existence of stable, localized solutions in the model The bound ()5.9p can be saturated for 
/3 = 1, if all the perfect squares in the second line in ()5.2j) vanish, leading to the Bogomol'nyi 
equations: 

„ , n(a — 1) „ , , 

fi+ ^ ^ V i = 0, (5.10a) 

(na-m) 

/2 + /2 = , (5.10b) 

P 

— + if? + f,-l) = 0. (5.10c) 
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The reduced action, S, can be obtained by replacing in the expression for the energy E in 



equations: 



+ Uq by u'^ 



^3 - ^0- 



Varying S gives the cyhndrically symmetric field 



{pa',y = 2as\f\'-2fi 



2/f(a-l) + 2/|(a- 



n 



-ipf[)' = h 



n 



P' 



P 



(5.11a) 
(5.11b) 
(5.11c) 
(5.11d) 



where ' = d/dp. These equations do not depend separately on uq and u^, but only on 



the Lorentz-invariant combination 



cUq, and therefore any solution determines 



actually a whole class, which is its Lorentz orbit corresponding to boosts ()4.3|) . Since we 
are interested in finite energy solutions of Eqs. ()5.1H) . in the following we shall assume that 
cu^ > (space-like or null classes in the terminology of the previous Section). 

We note that the change of variable a(p) — > nil — a(p)) achieves that the field equations 
do not depend separately on n and m, but only on the "helicity number" v = n — m. 
However, due to the boundary conditions, solutions of the field equations do not depend 
only on but separately on n, m. For this reason we prefer to use the parameterization of 
the equations explicitly containing both n and m. 

One can directly check that any solution of the first order Bogomol'nyi equations ()5.in|l 
also solves the field equations ()5.1H1 if /? = 1. In the next Section we shall briefiy review 
some relevant solutions of the Bogomol'nyi equations ()5.1()|1 . 



B. Local behaviour 

Eqs. (15.111) admit a 4-parameter family of local solutions regular at the origin, p = 0: 

a = a^'^p^ + 0(p2™+2) , as = af + Oip^^-'') , (5.12a) 

/i = + 0(p"+^) , h = ft^p"^ + 0(p-+2) , (5.12b) 

where a^'^\ af\ f'f'\ f^"* are free parameters. Depending on the relative magnitudes of the 
values of the parameters (3 and uj = \/uj^ with respect to -\/2, the various possible asymptotic 
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behaviours of the fields as p — oo for to > can be inferred from the following: 

„.l^^^e-v-._£;il^f^,..., ,5.13a) 

«3 = + + . . . , (5.13b) 

7p (1 - 2^2) p 

f2 = D^ + ... , (5.13d) 

VP 

where A, B, C, D are free parameters, and the resonant values (cj^ = 1/2, etc) have been ex- 
cluded. The dominant asymptotic behaviour of a, a^, fi is given by the smallest exponential 
in (jnSl). 

If ct;2 = then the Gauss-law-type equation for as ()5.11b ) (recall that ao = as) decouples 
from the field equations ()5.11b -d). The latter then become identical to the much studied 
"minimal" (i.e. static and untwisted) semilocal vortex equations. Let us recall some known 
facts about this "minimal" semilocal system ()5.11b -d). If /2 = 0, then Eqs. ()5.11b -c) reduce 
to the vortex equations in the N= 1 Abelian-Higgs model, in which case solutions are known 
to exist for any value of (3. The behaviour of these ANO vortices near the singular points 
is immediately obtained from Eqs. ()5.12|1 and ()5.13|) (with B = D = 0). In addition, 
Hindmarsh j?] and Gibbons et al. IiqI have shown the existence of n one-parameter families 
of solutions of Eqs. ()5.11b -d) also with /2 7^ 0, but only for the special value /? = 1. In 
fact, these solutions actually fulfill the first order Bogomol'nyi equations ()5.10|) . and they 
are sometimes referred to as "skyrmions". The asymptotic behaviour for large p of these 
solutions is radically different from that of the ANO vortices, whose fields approach their 
vacuum values exponentially, whereas the "skyrmions" exhibit an asymptotic power law 
behaviour. In fact for = the asymptotic behaviour of the fields in Eqs. ()5.1H) is given 
as: 

1 \ ^e-^P 



/II \ e-^P 

A = + (5.14) 
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where x = p/w with w, the width of the magnetic flux tube, being an arbitrary parameter. 
The magnetic field, falls off like 'uP'^ j p^^^'^ . This is quite contrary to what one could 
expect, namely that the flux be confined within a region of size ~ l/"\/2 determined by the 
length scale of the massive vector boson. It should be stressed at this point that although 
the asymptotic expansion (j5.14j) is valid for any value of /?, global solutions of (jS.lip for 
uj = satisfying the vortex boundary conditions ()5.12|) . ()5.14|) are known to exist only for 
(3 = 1. 

In this paper we shall give convincing numerical evidence that by relaxing the condition 
u = 0, n distinct one-parameter families of (twisted) solutions of Eqs. fl^.llj) exist not only 
for /3 = 1, but also for any value of /3 > 1. 



VI. UNTWISTED AND CHIRAL SEMILOCAL VORTICES 

Before considering new solutions, it is instructive to reproduce numerically some of the 
known solutions of the field equations (jS.llj) and to briefly review some of their properties. 
The numerical method we have used is a suitable version of the "multi-shooting" procedure 
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successfully applied to produce high precision data for the 't Hooft-Polyakov monopoles 
to which paper we refer for further details. 

The simplest solutions are the ANO vortices corresponding to = 0, 03 = 0, /2 = 0. 
We exhibit a sample ANO vortex profile function, a = Qano, fi = /ano; in Fig- 2, and also 
tabulate some numerical values of the relevant parameters of the ANO solutions in Table I. 
The energy of the ANO vortices increases monotonically with jS and diverges in the London 
limit /? — > 00. Solutions with (3=1 are distinguished by the fact that they also satisfy the 
first order Bogomol'nyi equations ()5.10|) . 

Next, as already discussed, for the special value (3=1 there are n one-parameter families 
of solutions with /2 7^ generalizing the corresponding ANO vortex. These solutions, 
dubbed "skyrmions", satisfy the Bogomol'nyi equations ()5.10|) . therefore they all have the 
same energy. As noted in [7[, the three Bogomol'nyi equations imply that 

/2 = ^/i ' (6-1) 



with X = p/w, so that ()5.10|) reduce to two equations for fi and a. For a given value of 
the winding number n, the helicity number v = n — m can assume n values, v = 1, 2, ... n. 
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1 1.5 

ln{l + p) 




FIG. 2: The profile functions for the ANO vortex with /? = 2, n = 1 (left) and for one of the 
"skyrmions" with n = l,m = 0, (/?=!) (right). 



which label n different solutions of the Bogomol'nyi equations for any given value of the 
"size parameter" w. The existence of these solutions was established in j^Q. We exhibit 
a typical solution in Fig. 2. 

In the limit w —>■ the "skyrmion" solution tends smoothly to the (3=1 ANO vortex. 
Since we shall make a connection with the large size limit of the "skyrmions", w oo, 
we recall some of the main points of the analysis of Ref. . To study this limit, it is 

convenient to introduce the variable 

X = ln(/) + lln(l + x-2^), (6.2) 

which satisfies the equation 

Now one can see that for large values of w the solution of Eq. ()6.3|) is well approximated by 

which tends to x = as oo, so as the width of the vortex tends to infinity the 

"skyrmion" solution approaches the following configuration more and more: 
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For the limiting configuration ()6.5p one has /f + /| = 1, which means that the amphtude 
of the Higgs field is fixed, <l>^$ = 1. Taking into account the U(l) "phase symmetry" one 
reahzes that the Higgs field $ takes its values in the target space of a CP^ nonlinear a- 
model. For this reason the "skyrmions" in the infinite width limit approximate CP^-lumps 
[3]. It is worth pointing out here that the non-trivial part of the limiting configuration ()(j.5|l 
depends only on z/ = n — m. For a fixed value of v there exist an infinite number of different 
"skyrmions" with different winding numbers n = z/, z/ + 1, . . . whose energies are different. 
Although solutions with n > u do not have the same boundary conditions at the axis as 
(lf).5|l . away from the origin all of them approach for w ^ oo the same configuration 

Finally one may wonder about the possible existence of such "skyrmions" with uj = also 
for /3 7^ 1. Just by counting the number of free parameters: there are four at p = in ()5.12|1 
and also four at p = cxo in ()5.14|) . one sees that one cannot hope to find a one parameter 
family of solutions similar to those in the (3 = 1 case. Although this parameter counting 
argument leaves open a possibility of finding some isolated "skyrmion" solutions, here we 
can only confirm the negative results of Ref. Q, since our numerical searches yielded no 
solutions with a; 7^ for /? 7^ 1. 



VII. TWISTED STRING SOLUTIONS 



In this Section we present some numerically constructed solutions of the cylindrically 
symmetric field equations ()5.11|) for (3 > 1. For a given value of the winding number 
n, the new solutions comprise n one parameter families labeled by the twist parameter, 
UJ. Exploiting the boost symmetry gives finally two-parameter families of twisted vortices 
labeled by u and uq. We also explore the phase space of these twisted vortices, and we believe 
that our mostly numerical results provide strong evidence in favour of their existence. In 
order to obtain numerical solutions by our shooting and matching techniques it is necessary 
to have a "reasonable" starting point in the 8-dimensional parameter space. In particular 
one has to avoid that the numerical procedure converges to the ANO solution, which is 
always present. A very natural idea is to search for new classes of solutions bifurcating with 
the ANO vortex. If such bifurcating solutions exist, they will provide a good starting point 
in the parameter space. 
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TABLE I: Numerical values of the shooting parameters f[^^ and a^^), 
the energy for the ANO vortex and the bifurcation parameter corre- 
sponding to the branching off of the twisted solutions. 

/3 ^b(/3) /f^ a(^^ E/27r 

1 0.500000 0.853178 1.000000 

2 0.329886 0.616573 1.099352 1.156761 

3 0.513051 0.696864 1.284290 1.260901 

4 0.661658 0.759706 1.438747 1.340595 

5 0.790405 0.811948 1.574136 1.405773 

6 0.905687 0.856955 1.696163 1.461215 

7 1.011028 0.896661 1.808164 1.509620 

8 1.108629 0.932290 1.912281 1.552676 

9 1.199978 0.964672 2.009991 1.591515 

10 1.286135 0.994398 2.102359 1.626935 



A. Bifurcation analysis 

The first step in the bifurcation analysis is to linearize Eqs. (jS.llj) around the ANO 
solution with a = aj^^o, fi = /ano cts = 0, /2 = 0. If now /2 ^ 1, then to first order in 
the amplitude /2 one has as = 0, since /2 = implies = 0, and /2 appears quadratically 
in Eq. (jS.llb ). Therefore 03 = 0(/|), and the linearized field equation ()5.11t i) for the /2 
amplitude can be written as the following (Schrodinger-type) equation: 
1 



- fLo) 



/2 = -^V2. (7.1) 



This equation corresponds precisely to the eigenvalue problem studied by Hindmarsh in his 
linearized stability analysis of the ANO vortex embedded into the semilocal theory Q. In 
agreement with his results, for any (3 > 1 we find a single normalizable solution, /2, with 
a negative eigenvalue, — cu^ = —uj'^{(3,n,m). In Table 1 we list for n = 1, m = the 
values of uj^iP) for some values of /3. Just as in Ref. Ql we do not find any non-trivial 
normalizable solutions of Eq. (j7.1|) when f3 < 1. In Ref. where 003 = 0, the existence 
of normalizable solutions (bound states) of Eq. ()7.1|) for /3 > 1 is the sign of the instability 
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of the embedded ANO vortex, whereas their absence for (3 < 1 indicates its hnearized 
stabihty. The instabihty of the embedded ANO vortex can be interpreted as the magnetic 
flux spreading out to infinity. 

In our case ^ 0, and the very same bound states indicate the existence of new branches 
of solutions bifurcating with the corresponding ANO vortex at uj = uJh{j3,n,m). The bifur- 
cation parameter ujy,{(3,n,m) is depicted in Fig. El in function of (3 for some of the smallest 
values of n and m. To each curve u!^{l3,n,m) in Fig. El there corresponds a new family of 
solutions bifurcating with the corresponding (i.e. of the same values of (3, n) ANO vortex 
for uj = uj\:){f3,n,m). Based on our numerical results there is little doubt that for any (in- 
teger) value of n a discrete family of curves uj^{l3, n, m) exist with m = 0, . . . ,n — 1. The 
new (static) twisted vortex solutions are therefore expected to exist for any value of (3 > 1, 
and for a given winding number n they comprise a two parameter set labeled by the twist 
parameter, u, varying continuously, and by the discrete index m = 0,1, . . .n — 1. Instead 
of the twist, u, each of the new solution branches can also be parameterized by /2™''. For 
m = 0, f!^^ is the value of the condensate at the axis of the twisted vortex. 

As can be seen in Fig. IHl all curves uj^{f3,n,m) pass through zero at /? = 1. This point 
corresponds to the bifurcation of the "skyrmions" in the limit of vanishing size parameter 
w with the corresponding ANO vortices. The corresponding zero mode solutions of the 
eigenvalueproblem describe deformations of the "skyrmions" with respect to w in the vicinity 
of 1(7 = [l]. Our numerics indicate that in the opposite limit, (3 ^ oo, one has ujh{(3, n, m) — *• 
oo. 

B. New Solutions 

As exposed in detail in the previous subsection, the known "magnetic spreading" insta- 
bility of the embedded ANO vortices in the EAH theory indicates in fact the bifurcation 
phenomenon, in which new (twisted) solution families branch off from the ANO family. In 
Subsection A we have already constructed new solutions in the linear approximation, but 
it remains to construct them in the full nonlinear theory. In view of the mathematical com- 
plexity of Eqs. (jS.llj) we do not have much hope to obtain analytical solutions, although 
an existence proof using functional techniques appears to be feasible. In this paper we 
present exclusively numerical results, which in our view provide convincing evidence as to 
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ujl{l3,n,m) 



n = 2, m = 



n = 1, m = 



n = 2, m = 1 
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FIG. 3: Left: The bifurcation parameter, uj^{[5,n,m), as a function of (5. Right: the value of the 
condensate as a function of the twist lo for twisted solutions with /3 = 2, 4, 9, and n = 1. 



/3 = 9 



/3 = 2 




FIG. 4: The rest frame energy E/2'k (left) and the rest frame current X/27r (right) plotted against 
the value of the condensate f!^^ at the origin for the solutions with n = 1 and /? = 2, 4, 9. 



obally regular twisted solutions. Adapting the multi shooting procedure 



the existence of g^ 

described in Ref. j2^ and tested for the ANO vortices (see Table I), we have succeeded to 
obtain numerically globally regular solutions of Eqs. ()5.11|) corresponding to twisted, current 
carrying vortices. To obtain the new solutions we start close to the bifurcation point, i.e. 
we choose values of the shooting parameters close to those of the bifurcating solution, and 
then we let vary the value of uj in the neighborhood of uj^{l3,n,m) to obtain "nontrivial" 
(i.e. different from the ANO vortex) solutions. 
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In fact for our solutions uj varies in a finite interval, < c<j < uj\^{f3,n,m) (for n > 2 the 
situation is actually slightly more complicated; see below). As oj tends to uj\,{(3,n,m), the 
functions as — > 0, /2 — > 0, while a ciano, /i /ano, i-e. the new class bifurcates with 
the corresponding ANO vortex. The limit — > is somewhat more complicated. Instead 
of approaching an untwisted vortex with = in this limit, solutions for u; — > approach 
pointwise the exterior part of the "skyrmions" in the large size limit {w — > oo), without ever 
attaining it. In other words, the point = does not belong to the phase space of twisted 
vortices. This limit will be analyzed in some detail in Subsection E. 
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FIG. 5: Twisted semilocal vortex solutions for n = 1, /3 = 2. 



Every solution of Eqs. ()5.1H) can be Lorentz-boosted, giving rise to a whole family of 
vortices, labeled by ujq, with a given u"^ = uj^ — uJq. Passing to the rest frame, where ujq = 0, 
we observe that the rest frame energy, E, decreases for decreasing values of u or equivalently 
for increasing values of the condensate f2^^ at the origin. This is illustrated in Fig. |5 where 
E/2tt is plotted in function of As the twist parameter decreases, the current and the 
condensate value at the center of the fundamental vortex increase; see Fig. El We have also 
plotted in Fig. E] the rest frame current, X/27r, (introduced in Eq. ()5.8|) ) in function of the 
value of the condensate at the origin, 

For solutions with n < 2 the decrease of the energy is monotone, so all such twisted 
vortices have lower energy than the corresponding ANO ones. For n > 2 the maximal value 
of the twist parameter can be different from the bifurcation value. For example for n = 3 we 
find that a;max(/5, 3, 0) > u;b(/5, 3,0) (this is recorded in Table IV), and so for a given value 
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of the twist parameter in the interval uj\^ < lo < cjmax there are two distinct twisted vortices. 
In this case /g'^'* can be chosen instead of on to label the solutions, since is always increases 
monotonically as one moves along the solution family. For sufficiently large values of /2°'* 
both uj and E start decreasing, and the energy E of the twisted vortex becomes smaller 
than i?ANo while u eventually tends to zero. The rest frame energy of the twisted vortices 
is therefore bounded as 

Eano{P, n) > E{(3, UJ, n, m) > 2?™ + u; |X|, (7.2) 

for n = 1,2, while for n > 3 the first inequality here is valid when u is sufficiently smaller 
than ciJb- The energy of the new solutions with = 1, 2 is thus always smaller than that of 
the ANO vortex. Therefore we expect the new solutions of winding number n = 1 with a 
nonzero current to be stable. 




lu(l + p) ln(l + p] 



FIG. 6: Twisted semilocal vortex solutions for (3 = 2 and n = 2, m = 0, w = 0.1 (left) and 
n = 2, m = 1, u) = 0.215 (right). 

In Figs. El we depict two sample twisted vortex function profiles for u = 0.3 and also for 
a rather small value of the twist, u = 0.017 (with (3 = 2 and n = 1). In Table II we list 
the relevant parameters of some of the twisted vortex solutions for n = 1: the twist, u, the 
values of f[^\ f^°\ a*^', af^ parameterizing the solutions at the origin, the rest frame energy, 
E, the invariant charge Q, and the rest frame value of the current, I = uQ. 

Next we turn to solutions with winding number n > 1. In this case one finds n distinct 
one-parameter solutions indexed by m = 0, 1, . . . , n — 1, or equivalently by the helicity index 
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TABLE II: Solution parameters for /3 = 2,4, 9 and n = 1. 



/3 


LO f^'^ a(2) af E/2TT Q/Itt I/2tt 


2 


0.3 1.010824 0.386212 0.546391 0.124789 1.155180 0.170624 0.051187 
0.2 0.767653 0.705426 0.363046 0.422267 1.131355 0.902637 0.180527 
0.1 0.534459 0.864141 0.204976 0.663061 1.083558 3.015082 0.301508 
0.017 0.239896 0.972041 0.051176 0.906725 1.020255 29.91849 0.508614 


4 


0.6 1.305601 0.428496 0.675653 0.107975 1.336898 0.096934 0.058161 
0.3 0.799348 0.837937 0.354284 0.496610 1.237417 0.869072 0.260722 
0.1 0.462641 0.947953 0.154324 0.754081 1.106005 4.140499 0.414050 
0.03 0.272278 0.981800 0.062910 0.889577 1.039866 18.60003 0.558001 


9 


1 1.677185 0.570981 0.808482 0.129739 1.573669 0.087067 0.087067 
0.6 1.113827 0.849213 0.515067 0.379056 1.443952 0.387808 0.232685 
0.3 0.734885 0.939568 0.300118 0.589641 1.273394 1.130746 0.339224 
0.1 0.435594 0.979182 0.136570 0.785287 1.116579 4.663623 0.466362 
0.08 0.394789 0.982893 0.116465 0.812686 1.097425 6.154640 0.492371 



TABLE III: Parameters of the fundamental solutions for [3 = 2 and n = 2, m = 

u /f a(2) af EI2'K Q/27r J/27r 

0.5 0.461753 0.743062 0.370423 0.500163 2.389129 0.459068 0.229534 
0.1 0.075237 0.990196 0.030212 0.951224 2.099932 4.633715 0.463371 
0.02 0.014362 0.999592 0.001576 0.996961 2.021527 26.21910 0.524382 



V = 1, . . . ,n determining the value of the angular momentum in Eq. ()5.6|) . In Figs. IHl we 
depict an n = 2, m = twisted vortex with u; = 0.1 and also another one with = 2, m = 1 
and ijj = 0.215. For m > both fi and /2 vanish at the origin, i.e. there is no condensate 
there. While the /? = 1 "skyrmions" with different values of m are all degenerate in energy, 
for twisted vortices with (3 > 1 this degeneracy is lifted. In fact, solutions with m > can 
be viewed as 'excitations' of the 'fundamental' solution with m = 0, since they have higher 
energy, this is illustrated in Fig. [7| 
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FIG. 7: Rest frame energy E/{2'k) versus the rest frame current T/(27r) for two solutions with 
n = 2, m = 0, 1 for (3 = 2. For m = and for large currents one has E/ (27r) — > 2. 

TABLE IV: Parameters of the fundamental solutions for /? = 2 and n = 3, m = 

w /f^ a(2) af^ E/2iT Q/2it T/2tt 

0.71 0.389410 0.122714 0.647036 0.018799 3.716103 0.010221 0.007257 
0.73 0.353581 0.380272 0.587955 0.167458 3.717523 0.094779 0.069189 
0.4 0.064950 0.965118 0.101566 0.878371 3.459984 1.215289 0.486116 
0.13 0.009128 0.998709 0.007153 0.988680 3.167331 4.627920 0.601630 



We tabulate a few values of the pertinent parameters of some fundamental twisted solu- 
tions with n = 2, 3 in Tables III and IV, and of some excited ones with n = 2 and m = 1 
for /3 = 2 in Table V. 

In the hmit u cfb(/3, n, m) solutions with different values of m tend to the same ANO 

TABLE V: Parameters of the excited solutions for /3 = 2 and n = 2, m = 1 

/f^ f2^ a(^) af^ E/2iT Q/2it T/2tt 
0.2 0.559894 0.721888 0.649659 0.374412 2.403820 1.288194 0.257639 
0.1 0.343013 0.976403 0.641370 0.715442 2.321518 5.587675 0.558768 
0.04 0.212548 1.054153 0.629574 0.873387 2.241088 20.31268 0.812507 
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vortex, whose energy does not depend on m. In the opposite hmit, for small values of the 
twist and for large currents, the energy of the fundamental solutions approaches the universal 
lower bound E = 2Tm irrespectively of the value of (3, while for the excited solutions this is 
not the case, since their energy grows with (3 and diverges in the London limit (3 ^ oo for 
any I (for a more detailed discussion of this limit see the next Subsection). 

We now describe some interesting features of the stationary case. As has already been 
discussed stationary solutions can be obtained from the static ones by applying Lorentz 
boosts ()4.4|) . For static solutions the global charge, momentum and angular momentum 
densities are all vanishing, Jq = = = 0. Applying a Lorentz boost to a static 
solution generates a stationary one whose fields are obtained by simply replacing in Eq. 1)5.111 
oj^ = uj and cuq = by cuq = — co'sinh(7), oj^ = cosh(7)co'. These stationary solutions have, 
however, some different physical properties, for they have non- vanishing momentum, angular 
momentum, and global charge. It is instructive to see how this comes about explicitly. 

Let us first remark that the energy momentum tensor of the ANO solutions is boost 
invariant, since its only non-zero components are T^q = T^^, T^p, T'^^. When we compute 
T'^^, for the static twisted vortices, we find that = T^^ = T^^ = 0, but otherwise the 
remaining 7 components are all non-vanishing and describe a complex system of internal 
stresses and pressures in the vortex. The boost invariance of this system is manifestly 
broken by non-zero values of T*}, — T^^ and T^^, which can be represented in the form 

n-n = |ji+---, n = ^^'i+---, (7.3) 

where dots stand for total derivative terms that vanish upon integration. If we now apply 
a Lorentz boost ()4.4)1 . then in view of the tensorial transformation law, T^iy — > T^^, = 
{dx" /dx^){dx^ /dx'')Tai3, the densities of the energy, momentum and angular momentum 
become 

n = n+sinh2(7)(n-n), 

f° = -sinh(7)cosh(7)(n-T^J, 

= sinh(7)r; . (7.4) 

Integrating Eq. ()7.4|1 using Eq. ()7.Hj) and / j| (fx = 4X3 = AuQ one can express the effect 
of a boost on a static solution as E E = E + 2{ijJq)'^Q, P = 2ujquj^Q, J = —2ujqI'Q, 
which completely agrees with Eqs. ()5.2|) . ()5.6p as it should. Similarly the charge and current 
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densities generated by such a boost are j'q = — sinh(7)j| and j| = cosh(7)j|, which gives 
upon integration Eq. ()5.5|) . 
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FIG. 8: Left: The energy and current densities, Tq and j'l, for a static solution with (3 = 2, 
uj = 0.19, and n = 1. Right: Energy densities, Tq , of the corresponding stationary solutions for 
some values of the boost parameter, 7. 

Summarizing, twisted vortices are not Lorentz boost invariant. The densities of the 
energy, momentum, angular momentum, global charge and the current of a boosted static 
solution are given by linear combinations of T\ and j|, up to terms that vanish upon 
integration (see FiglH)). These stationary solutions can be labelled by the value of the boost 
parameter, 7, or equivalently by ujq, or equivalently by their momentum P. 

C. The (3^ 00 limit 

It is instructive to discuss in some detail the (3 = 00 limiting theory, in which case the 
scalar fields are constrained as 

<f^<^> = l/iP + l/2p = l, (7.5) 



when the SU(2) EAH model reduces to the gauged CP^-model. It is convenient to param- 
eterize the scalar fields as: 

fi = cos9, f2 = sm9, (7.6) 
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and then the field equations can be written as: 



-{pas')' = 2[a3 - cos^ 6] , 



/ a' 
p{ - 
\P 

1 



P 



ipO' 



cos 



n 



w'(2a3 - 1) 



nu 



P' 



2a 



n + m 



n 



sin(2e) 



(7.7a) 
(7.7b) 
(7.7c) 



The local behaviour of the fields near p = is not analytic in p, it is given by 



a = (a*-^-* In p)p^ + . . . , as = a?^ + . . . , 

n 



e = cp 



yjn'^—m'^ 



+ . . 



(7.8) 



For large values of p the asymptotic series of the function 9 is easily seen to be 



-ujp 



e = n/2-D + ... 

Vp 



(7.9) 



while the various possible asymptotic behaviours of a, 03 for p ^ 00 depending onu ^ 1/ -\/2, 
are still given by the formulae ()5.13b .b). The profile functions of the current carrying vortex 
solution of Eqs. ()7.7|1 are qualitatively similar to those of Eqs. (15. Ill) . for a sample profiles 
see Figs. El In Table VI we give the relevant parameters of some solutions with winding 
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FIG. 9: Semi-local vortex solutions for n = 1, /? = 00. 



number n = 1 for various values of the twist, u. For any winding number n the energy 
of the fundamental {m = 0) (3 = 00 solutions remains bounded for any finite value of uj. 
We emphasize that our numerical evidence does suggest ujy,{oo) = 00. i.e. the parameter uj 
varies in the interval (0, 00). 
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TABLE VI: Parameters of the n = 1, /3 = oo solutions 

oj f['^ a(2) of E/2TT Q/27r J/27r 
5 6.030722 1.812331 0.064034 2.991526 0.018523 0.092613 
1 1.442394 0.640988 0.366642 1.753574 0.274199 0.274199 
0.1 0.417767 0.125292 0.804925 1.124363 5.046921 0.504692 
0.01 0.155052 0.022236 0.958536 1.017715 78.02346 0.780235 



On the other hand, the energy of the excited (m > 0) /? = oo solutions diverges, since 
the (m — na)^/|/p term in the energy density ()5.2j) causes a logarithmic divergence in the 
total energy as p — 0. This behaviour is very similar to the case of the f3 ^ oo limiting 
ANO vortex, whose energy also diverges logarithmically. 

D. Qualitative description of the phase space 

The phase space of the fundamental (m = 0) twisted vortex solutions with winding 
number n = 1 is depicted in Fig. ^[ using the explicit parameters appearing in equations 
(I5.11|l . jS, uj, and the magnitude of the condensate at the origin, f2^\ as coordinates. In 
Fig. El a one parameter family of twisted vortices for a fixed value of (3 corresponds to a 
curve in the (/2°\ on) plane. The ANO vortex for some fixed value of (3 corresponds to a 
f3 =const. straight line in the /g^'' = plane. The curve t^b(/5) is also indicated in the /g'^'* = 
plane, where the bifurcation of the new solutions with the ANO vortices takes place. Three 
families of twisted solutions have been indicated as solid curves for the values jS = 2,4,9, 
together with the "skyrmion" solutions satisfying the Bogomol'nyi equations {(3 = 1). The 
curves terminate at the smallest values of u we could compute. They should all approach 
indefinitely the {ui = , f^"* = 1) line, without ever reaching it. 

We have also depicted the "energy landscape" of the fundamental solutions with winding 
number n = 1 in Fig. ^2 In Fig- a family of solutions for a fixed value of (3 corresponds 
to a solid curve in the {f2^\ uj) plane, and we have also indicated the corresponding ANO 
vortex, located at the endpoint of the curve in the /j^'' = plane where the new family 
bifurcates with it. The family of solutions for /5 = oo is represented by the solid curve in the 
f2^^ = 1 plane. It is clearly seen that as the value of u decreases from u}h{(3) towards a; = 0, 
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FIG. 10: Phase space of the fundamental n = 1 twisted vortices. 

the energy of all twisted vortices decreases towards the Bogomol'nyi bound, E = 2tt, without 
ever attaining it, irrespectively of the value of (3. An important point is that solutions with 
nonzero values of uj do not have a smooth oo ^ limit, i.e. the point (/g^'' = 1 ,uj = 0) 
corresponding to the approximate CP^ lump does not belong to the phase space as is 
indicated on Fig. [TT]by a circle. 

E. The zero twist limit 

As we have already remarked, the limit of zero twist, u; — is rather nontrivial. After 
a brief description of the behaviour of the solutions for small values of uj, we shall give a 
heuristic explanation of the observed behaviour without any attempt of mathematical rigor. 
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FIG. 11: Energy of the fundamental twisted vortices. 

One observes that when the twist parameter uo decreases with fixed /5, n, m, the solutions 
behave differently in a region containing the origin (interior region) and in the corresponding 
outside region. In the outside region the asymptotic behaviour of the solutions is dominated 
by the /2 amphtude due to its slowest falloff, but what is more, the solutions approach 
more and more the "skyrmion" configurations described in Section VI. At the same time 
the value of the current is getting larger and larger, in fact for n = 1 we expect that it will 
grow without bounds. We depicted the current in function of the twist parameter for n = 1 
in Fig. ^1 Irrespectively of the value of (3 the energy of the fundamental solutions (m = 0) 
approaches a universal lower bound which is nothing but the Bogomol'nyi lower bound, i.e. 
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FIG. 12: Left : the current, I/2-it, as a function of uj for n = 1, /3 = 2,4, 9, oo. Right : The energy, 
E/(2Trn), against the current for (3 = 2 and n = 1, 2, 3. 

for a; ^ we observe that X is growing without any apparent limit and 

E^27m, /f^l. (7.10) 

The excited solutions with m 7^ also approach the same limiting "skyrmion" configurations 
depending only on v. Unfortunately, with our numerical procedure we could not obtain 
solutions for values of uj smaller than ^ 10~^. Nevertheless these numerical results enable 
us to give a heuristic description of the — limit. It turns out that the solutions behave 
in this limit essentially in the same way as the large size "skyrmions". Let us introduce a 
new function C(p) as 

f2{p) = x-^C{p)h{p), (7.11) 

where C(0) = 1, x = p/w, and w is the "size" parameter. Similarly to the skyrmion case, 
we will now study the limit w 00. With 




(7.12) 
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the field equations ()5.1ip read 

1 



s - {xa'.y = 2as{ft + - V2 , (7.13a) 

sx(-] =2fl{a-l) + 2fl{a--), (7.13b) 
\x J n 

[k^f[y - ^"^^-^h} = h - /?(1 - fl - fl)] , (7.13c) 



// , l-2z/\ , 2z/ /// nil- a) 



UJ 



2 



C" + 2^ + \c'- — (^- ^ ^ C = — (1 - 2as)C, (7.13d) 

V/i ^ J X \fi ^ ) ) s 

where now ' denotes the derivative with respect to x. Assuming that 

hmcj2 = 0, (7.14) 

and that the derivatives remain bounded in an interval, T>, equations ()7.13b -c) imply on V 
for s — > 

= {fl + fl)a3-fl, = flia-l) + fl[a-'^), = /^(A^ + /| - 1)A, (7.15) 

whose solution can be parameterized with a single function B(x): 

/i = sin(e), /2 = cos(e), a3 = cos2(e), n{l - a) = u cos'^ (O) . (7.16) 

The function C{x) can now be expressed as 

C = x'^^ = x'^cot(e). (7.17) 
fl 



Inserting ()7.17|) to the equation for C{x) ()7.13t i) one obtains 

X A X 



e" + - e' = - (A + -) sin(4e) , (7.18) 



where A = lims_^o t^^^Z-s- Since the boundary conditions, 9(x) 7i/2 for x oo and 
Q{x) — ^ as a; 0, ensure that the fields /i,/2,a,a3 have the qualitatively correct global 
asymptotic behavior not only on V but by extension on the whole < x < oo, we seek 
solutions of ()7.18|) subject to these boundary conditions. There are now two possible cases: 
either A is unbounded or it is finite. In fact A must be bounded, since otherwise 0(x) could 
not satisfy the boundary conditions. A simple positivity argument shows that, unless A = 0, 
equation ()7.18|) does not admit solutions satisfying the boundary conditions. When A = 0, 
the solution is given by 
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We therefore conclude that in the = w — > oo hmit one has Mm.s^QUj'^/s = 0, and the 

solutions approach the following configuration 

VI + x^'' V 1 + x^'^ 1 + x^"" 1 + x^'" 

and furthermore C{x) = 1. This limiting configuration is indeed the same as the one 
obtained in the large size skyrmion limit ()6.5|) . What is interesting, is that this limiting 
behavior is universal, i.e. it does not depend on (3. Similarly to the skyrmions, the limiting 
configurations depend only on u, and so solutions with n = z/, i/ + 1, . . . have the same limit. 
Only the fundamental m = solutions approach the limiting skyrmion globally, since they 
have the same boundary conditions at the origin as ()7.20|) . and so their energy approaches 
in the limit the Bogomol'nyi bound, i.e. E = 2?™. The excited solutions (m 7^ 0) have 
different boundary conditions at x = 0, and they approach the limiting configuration in a 
domain not containing the origin. 

Finally we present a simple virial argument which shows that at least in the j3 = 00 
theory there are no nontrivial "skyrmion" solutions. From the energy functional ()5.2|) one 
can obtain the following virial relation for static configurations: 

I dp^- = j^ pdp\2uj'[flal + fl{l-a,Y]+(3{l-\f\^YY (7.21) 

which immediately implies that for = the (3 = 00 theory, when the the potential energy 
term cx /3 is absent, does not admit nontrivial finite energy solutions. Although for finite 
values of j3, the above simple virial argument does not exclude the possible existence of 
"skyrmions" for = 0, the analysis of the — limit has indicated that in fact irrespec- 
tively of the value of /5, /i + /| — > 1, therefore it is unlikely that nontrivial "skyrmions" 
would exist apart from the known ones for /3 = 1. 

VIII. SUMMARY AND CONCLUDING REMARKS 

In summary, we have presented a new class of straight, twisted string solutions of the 
SU(2) symmetric semilocal Abelian Higgs model. These twisted strings have finite charge, 
energy, momentum, and angular momentum per unit length, and they also carry a global 
current. The solutions are not Lorentz boost invariant in the t, z plane, and there exists 
a rest frame where their charge, momentum, and angular momentum vanish. The current 
does not vanish, however, and so it is an essential physical parameter. 
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The twisted solutions exist for /3 > 1 and they can be labeled by the frequency 
and twist uj^, by the winding number n = 1,2,..., and by an additional integer, m = 
0,1,... n — 1. Related to these, there are four independent physical parameters of the 
vortex: the momentum P ~ ujq, current X ~ u^, magnetic flux ~ n, and angular momentum 
J ~ {n — m). For X — > the twisted solutions bifurcate with the embedded ANO vortices. 
The current seems to range in the infinite interval, at least for n = 1, and for |X| oo the 
solutions approach the CP^ lump(s) just as the "skyrmions" do in the large size limit. The 
energy in the rest frame decreases as |X| grows and, if m = 0, it tends to the absolute lower 
bound 27171 in the zero twist limit. The fundamental solutions with m = have the lowest 
energy for a given n and X 7^ 0. They have finite energy also in the London limit (3 ^ 00, 
whereas the energy of the excited strings with m 7^ diverges. As u 0, the solutions 
approach the limiting "skyrmion" configurations. 

Perhaps the most important property of the twisted strings is the fact that their energy is 
lower as compared to that of the ANO vortices. They are therefore energetically more favor- 
able to be produced during phase transitions, and so one would expect that the formation 
of such twisted strings is an important issue. This can be called 'spontaneous superconduc- 
tivity' and this may imply interesting applications in cosmological models. The new twisted 
vortices could also play a role in solid state physics in models of various superconductors 
with several order parameters (e-g- two-band superconductivity models containing a vector 
and two complex scalar fields j2ll23|)- In superfiuid ^He twisted vortices qualitatively similar 
to ours (Salomaa-Volovik vortices) have been experimentally observed P|. From the point of 
view of high energy physics the twisted semilocal strings should have their twisted Z-strings 
counterparts in the bosonic sector of Standard Model of Electroweak interactions. 

An important property of the superconducting semilocal strings it that, similarly to 
the topological ANO vortices, they cannot terminate. The embedded ANO strings do not 
have this property, since they can terminate on magnetic monopoles [^. Superconducting 
semilocal strings have a global current flowing through them, whose conservation implies 
that they should either be infinitely long or form closed loops. Such loops will in general 
be charged and will also have a momentum circulating along them, and so that there is 
a possibility that they could form quasi-stationary objects of vorton type. In fact there 
have been attempts to construct such closed flux loop solutions in the SU(2) EAH model, 
although the results obtained remain somewhat controversial for the time being since the 
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2J], has not been confirmed in the analysis of 



existence of such solutions, advocated in Ref. 
Refs. 0- 

Maybe the most important open problem related to the new solutions is their stability. 
Since for a fixed value of the current the fundamental twisted solution have the least energy, 
we expect them to be stable in the semilocal model. However, a detailed stability analysis 
is necessary to make definite conclusions. 

Another interesting issue is to investigate the interactions between twisted vortices. While 
type II ANO vortices (with /? > 1) always repel each other, for our twisted vortices the 
situation seems to be more complicated. Since we have superimposed vortex solutions for 
n> 1, we can ask if it is energetically favorable for them to break up into widely separated 
constituents with n = 1. The curves in FiglT^ where the ratio E/{27m) is plotted against 
X, indicate that the interaction is repulsive up to a certain value of X, however, at least for 
n = 2 it becomes attractive for large enough currents. For the n = 2 vortex a clear crossover 
takes place around X/27r ^ 0.51. This suggests that the n = 2 superimposed vortex for large 
enough currents becomes stable with respect to breakup into n = 1 vortices, unlike the type 
II ANO vortices. If the interaction of separated vortices has also an attractive phase this 
could have important physical consequences. Clearly a more detailed analysis of this issue is 
necessary, which is, however beyond the scope of our paper. It is also worth mentioning that 
in a different two-component Ginzburg-Landau model a similar non-monotonic behaviour 
of E/n has been reported in [2^. 

Before finishing our discussion, we would like to comment on the fact that the vortex 
charge Xq, momentum P, and angular momentum J vanish in the rest frame, but the cur- 
rent X3 does not. This can certainly be explained using the tensorial lows of transformation 
of the corresponding densities, as was done above. At the same time, it might be worth 
illustrating the same fact also in the framework of a simple kinematic model, viewing the 
superconducting string as consisting of positive and negative elementary charges of the same 
mass in stationary motion. The positive charges spiral around the vortex axis clockwise in 
the positive z direction, while the negative charges spiral with the same speed but counter- 
clockwise and in the negative z direction. The positive charges have thus positive momentum 
and angular momentum, while the negative charges have negative momentum and angular 
momentum, but the current is always oriented clockwise in the positive z direction. In the 
vortex rest frame the volume densities of the positive and negative charges are thus identical. 
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implying that the total charge density, momentum, and angular momentum vanish, while 
the total current does not. If the vortex is boosted in the positive z direction, say, then, 
due to the Lorentz contraction, the volume density of positive charges will become higher 
as compared to the density of negative charges. As a result, the total charge, momentum 
and angular momentum will no longer be zero. 
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